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Abstract 



CsJ ' We extend a procedure for construction of the photon position operators with 

transverse eigenvectors and commuting components [Phys. Rev. A 59, 954 



(1999)] to body rotations described by three Euler angles. The axial angle 
can be made a function of the two polar angles, and different choices of the 



o . 

' functional dependence are analogous to different gauges of a magnetic field. 

Symmetries broken by a choice of gauge are re-established by transformations 
within the gauge group. The approach allows several previous proposals to be 



related. Because of the coupling of the photon momentum and spin, our posi- 



^ ■ tion operator, like that proposed by Pryce, is a matrix that does not commute 



with the spin operator. Unlike the Pryce operator, however, our operator has 

commuting components, but the commutators of these components with the 

total angular momentum require an extra term to rotate the matrices for 

each vector component around the momentum direction. Several proofs of 

the nonexistence of a photon position operator with commuting components 

are based on overly restrictive premises that do not apply here. 
PACS number(s): 03.65.Ta, 14.70.Bh, 42.50.-p, 03.65.Vf 
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I. INTRODUCTION 



Since the early days of quantum mechanics, it has been claimed that there is no pho- 
ton position operator with commuting components, and hence that a basis of its localized 
eigenvectors, |r'), does not exist As a consequence, it is widely held that there is no 
coordinate-space photon wave function, ip (r') = (r'\tp) [0. Wave functions for photons in 
momentum space are commonly used, and it is there where the position operator would be 
defined. While it may generally be possible to use a second-quantized formalism in which 
the fields are operators in Fock space [@,[|, the usefulness of a coordinate-space photon wave 
function for describing such phenomena as photon interference is sufficiently well recognized 
that a number of authors have introduced versions of effective spatial wave functions |||| . 
The existence of a position operator may also be important for the consistency of some work 
in the second-quantized formalism. 

A number of authors 0-0] have discussed and to some extent resolved well-known prob- 
lems of Lorentz covariance and causality that exist for any position operator in relativistic 
quantum theory. However, there are additional problems for massless particles of spin one or 
greater, and a number of "proofs" have been presented of the nonexistence of a photon po- 
sition operator in the literature |T3|-|I8f . In spite of these, one of us has recently produced a 
counter example: a Hermitian position operator with commuting components and localized 
transverse eigenvectors [19]. However, the asymmetric, singular nature of the operator was 
puzzling, and questions concerning its compatibility with work concluding the nonexistence 
of such operators were not fully addressed. 

In the present paper we attempt to resolve these issues. Our principal tool is a gener- 
alization of the new position operator to include an arbitrary axial rotation, which may be 
a function of the other two Euler angles that parameterize the body rotation. This gener- 
alization provides insight into the geometry underlying the operator and allows us to unify 
several previous approaches. 



While position operators can be defined in several ways, the one in Ref. [IS], henceforth 
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to be referred to as I, was constructed by requiring its components to have eigenvectors 
transverse to the momentum and in the directions of the momentum-space polar unit vectors. 
The position operator thus obtained takes the form of a 3 x 3 matrix and will be referred 
to here as r^°\ where the underscore denotes the matrix character and (0) refers to use 
of the spherical polar basis vectors 8 and 0. The momentum-space operator r*- -* can be 
expressed by a transformation that rotates the components of the photon state function in 
momentum space to a fixed photon reference frame, differentiates, and then rotates back to 
the lab frame. Since two angles suffice to specify the direction of the photon momentum p, 
the rotations require only two independent parameters. 

Here we follow a more general approach by including, in addition to the polar angles 9 
and 0, the axial Euler angle, x- Let the unit vectors along the Cartesian axes be ei, e 2 and 
e 3 . A rotation about p by an angle x permits an arbitrary choice of the transverse unit 
vectors e p i and e P 2 , which are obtained by a rotation from the fixed unit vectors ei and e 2 
in the given reference frame. The same rotation is designed to take e 3 to the momentum 
direction p = p/p, but p depends only on the polar angles 8 and 0; it is independent of x- 
We show below that if x is independent of 9 and 0, the position operator is just as found 
previously However, we are also free to choose x to be a function x = Xp 4>) °f the 



polar angles. The momentum-space position operator that corresponds to such a choice is 

r = r(°) + Vx P (1) 

in units with H = 1. (It is understood that the second term on the rhs is multiplied by a 
unit 3x3 matrix. To simplify notation, we do not indicate unit matrices explicitly.) A 
change in the momentum-space gradient Vx P (#,0) is analogous to a gauge transformation 
of r. The singularity previously found in is analogous to that of a monopole string and 
commonly arises in electromagnetic-like gauge potentials. For any choice of Xvi L h as less 
symmetry than one would expect and under parity inversion it does not simply change its 
sign. However, as we show below, the symmetry is restored when one includes the full group 
of possible gauge transformations rather than only a single choice of gauge. 



In the case of massive particles of spin S, any component of the spin can be used to 
define an eigenstate basis. The basis functions corresponding to a space-fixed quantization 
direction span a rotationally invariant linear manifold and carry irreducible representations 
of the Poincare group. Because the momentum-space basis states for massive particles 
have directions in space that are invariant under infinitesimal translations, the position 
operator can be identified with the Hermitian generator iV of such translations, where V is 
the momentum-space gradient. Since zV is a differential operator proportional to the unit 
matrix, it commutes with spin matrices that represent space-fixed components of the spin. 
For massless particles, on the other hand, representations of the Poincare group are reduced 
to the two irreducible representations carried by helicity states in which the spin component 



along the momentum p is ±S ||20|| . Since the two helicity subspaces are invariant under the 
actions of the Poincare group, the allowed spin states for a massless particle with S > 1 form 
rotationally invariant manifolds with fewer than the 2S + 1 independent elements required 
for states quantized along a given spatial direction. As elaborated below (see Section VI), 
this is important in understanding the relevance of the nonexistence result of Newton and 



Wigner [TB|]. The basis vectors corresponding to the two allowed helicities are transverse to 
p and therefore require a modified position operator, one that is no longer independent of 
the spin. 

The generator of rotations is, as usual, 

J = -zpxV + S, (2) 

where the vector components of S are the spin-one matrices Sj whose explicit elements 
depend of the basis used for the state vectors (see discussion in the following section). The 
term — ip x V = L is required in order to generate a rotation of the argument of the wave 
function, whereas S generates the rotation of its components. Our use of the sum (0) for 



photons is consistent with Bargmann and Wigner's |21[ proof that J is the generator of 
rotations for a particle of arbitrary mass. 
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For the position operator r to transform as a vector, it is traditionally required that 



J j i Tk 



iejkin , (3) 



where repeated indices are summed over. For massive particles, this commutator follows 
directly from the assumptions that components of the operator r commute with each other 
and with the spin matrices, and that they satisfy the canonical commutator relations with 
the momentum, 

\pj,r k ] = -i6 jk , (4) 
since then the components of either L or J rotate one component of r into another: 

Jji r k = [Lj,fk] = £jrnn r in[Pm r k] (5) 
i& ' jkmTm • (6) 

Whereas this clearly applies to the position operator zV usually chosen for massive particles, 
the photon position operators proposed by Pryce and one of us JOJ are matrices that do 



not commute with the spin. As we discuss in more detail in the following section, components 
of the Pryce position operator satisfy the commutation relations (|3|) but do not commute 
with each other. Consequently, their eigenstates cannot have eigenvalues for more than one 
component of the position operator and thus cannot fully define the position of a photon. 

On the other hand, our position operator has commuting components, but at the cost of 
the commutation relation (^). For a given functional form of x P , an infinitesimal rotation of 
r results in an incremental axial rotation of the matrix rj for each component by d\ — dx P , 
and this yields the extra term derived below [see Eq. (|67Q1. Recent proofs of the nonexistence 
of a photon position operator and localized states have assumed components rj that satisfy 
Eq.([|) JT7| , |T8|1 , thus excluding the position operators discussed here. 

As in I, the position operator is constructed with the requirement that its components 
have transverse eigenvectors. It is related to the position operator for massive particles, 
whose components have eigenvectors with fixed directions in momentum space, by a trans- 
formation with the form of a spin rotation and a dilation. As discussed more fully in the 



following section, this similarity-like transformation includes a spin rotation D (0, 9, x) from 
fixed directions in the photon reference frame to the longitudinal and transverse directions 
in the lab. Because the longitudinal direction is independent of the axial angle Xi a family of 
transformations with different functional dependencies x = Xp ($> 0) are possible. We find 
that this family of transformations includes ones that relate the Shirokov and Lomont and 
Moses irreducible helicity representations to the Foldy form of the Poincare algebra. The 
use of Euler angles thus gives a unified approach to the relationships among these represen- 
tations and to the position operators with transverse eigenvectors obtained in the present 
paper. 

In Section II we briefly review the Foldy and helicity representations of the Poincare 
operators. We then describe the relationship of the angular momentum and boost operators 
to the Pryce |22| position operator. In Section III the results of paper I are reviewed and 
then generalized to an arbitrary functional dependence x = Xp 0) °f the axial Euler 
angle on the polar angles of p. The Poincare group, the relationship of the position operator 
to the angular momentum and to the boost operators, and the localized states are then 
examined in Section IV. The connection of our work to Berry's phase is treated in Section 

V, and the relation of string singularities to nonintegrable angles is stressed. In Section 

VI, the consistency of our position operator with nonexistence proofs is addressed, and its 
implications for the photon wave function are briefly discussed. Details and alternative 
derivations of some of the geometric arguments are presented in the Appendices. 

II. FOLDY AND HELICITY REPRESENTATIONS OF THE POINCARE 

OPERATORS 

Foldy [^3|] wrote down the form of the ten generators of inhomogeneous Lorentz trans- 
formations (the Poincare operators) for particles with spin S and mass m in the standard 
space-fixed representation. His operators also exist in the m —> limit, the only case con- 
sidered here, and in momentum space they comprise the momentum p, the Hamiltonian 
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H = pc, the total angular momentum operator J [Eq.@], and the boost operator 

K = ipV + p x S. (7) 



(Foldy's symmetrized form of the ipS7 term was criticized by Chakrabarti |[24|| . The correct 
form is related to the momentum-space normalization weight used for scalar products, which 
is discussed toward the end of this section. Note that the boost generator K is — N in 
Chakrabarti's paper ||24|| .) 



The standard helicity representation was introduced by Lomont and Moses [p5|] to provide 
a realization of the generators for the zero-mass case. The representation is distinct from 
the zero-mass limit of the Foldy representation, given above, because the carrier space has 
been split into invariant subspaces labeled by the two helicity components k = ±S. Vector 
components of the spin S are referred to the photon reference frame in which the momentum 



direction p is 63. Chakrabarti [24| showed that the zero- mass limit of the Foldy representation 



is related to the Lomont and Moses helicity representation by the unitary transformation 

F = UO LM U-\ (8) 

where Op is an operator in the Foldy representation, Olm is the corresponding operator 
in the helicity representation of Lomont and Moses, and the unitary transformation has a 
matrix representation U_ that we recognize as the spin rotation 

U = exp{-i64>-S). (9) 

The transformation @ rotates S3 directly into S • p. As above, 9 and are the usual 
spherical polar angles of p in momentum space and <f> is the unit tangent vector in the 
direction of increasing <fi. It is perpendicular to both e 3 and p . In the transformation 
(§) and its inverse, it is important to recognize that the unitary operator Eq.(|9|) does not 
generally commute with the momentum-space gradient V because of the dependence of 9 
and <f> on p. The Hermitian generators of the infinitesimal transformations of the Lomont 
and Moses helicity representation are thus given by the Foldy representation and the inverse 



of (D to be d|] p = J2jPj e j> H = pc, 



— LM 



— ip x V+ 



P + e 3 \ o 

i + cos^y — 



(10) 



and 



P x e 3 

1 + cosfl / — 



5, 



(11) 



where the matrix 5*3 can be replaced by k when restricted to the invariant subspace of 
definite helicity. 



A different helicity representation was obtained earlier by Shirokov |26j . After conversion 
to polar coordinates, his angular-momentum and boost operators are 



J Sh = -ip x V+ (d cot 9 + pj S3 



(12) 



and 



K Sh = ipV + 4> cot 6 S 3 . 



(13) 



The Bialynicki-Birulas ||27| , |28|| derived results equivalent to Eqs. (|l2l) and (^) by considering 



the Maxwellian momentum and angular-momentum tensors. They also found j28[ a unitary 
matrix W\ that relates operators in the Foldy and Shirokov representations: 



(14) 



We can express their W\ as the product of a transformation T from a Cartesian to an 
angular-momentum basis and another rotation 

D (0, 9, 0) = exp (-iS^) exp (-iS^O) (15) 

that takes e 3 into p : W\= D_ ((f), 0,0)T_ = T_D_^ (0, 6, 0) . The unitary transformation 

d i\ 



T 



V2 



i -i 

yo V2 j 
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(16) 



is required only because different bases for the momentum-space vectors have been assumed 
in the two representations. Specifically in the transformation (|I3D, the Foldy operator Of 
uses a Cartesian basis, in which the unit vectors e 1; e 2 , e 3 are represented by column vectors 
with elements (&,-). = 5jk, whereas Og^ uses an angular-momentum basis in which the 
complex unit vectors e re , k = ±1, 0, are column matrices with elements (e K ) k = S K ^-2 ■ They 
are given in the Cartesian basis by the columns of T (|16|) . The spin matrices are different in 
the two bases. In the Cartesian basis, the matrices Sj have elements (Sj) kl = —iejki whereas 
in the spin-1 angular- momentum basis the corresponding matrices are ( S3 ) = (2 — k) Ski, 
and is^p ± iS^ ) = V28k,i±i ■ The body rotation matrices, expressed in terms of the Euler 
angles 0, 9, % by 

D((f>, 9, x) = exp (-iSa0) exp (-iSrf) exp (-iftx) > ( 17 ) 

have correspondingly different elements in the two bases. Explicit values in the Cartesian 
basis are given in the following section and in Appendix A. Elements of D} 1 ' = T^ l D_T_ 
in the angular-momentum basis have the standard form D^, = exp (— ik§ — ik'x) d^, (9), 
where the real matrix dS 1 ' is [29| 

( cos9+ 1 v^sin^ cos9 - 1 ^ 

-v^sin^ 2cos0 -v^sinfl . (18) 
, cos 9 — 1 \pl sin 9 cos 9 + 1 , 



d (1) (0) = \ 



We will generally assume that the different representations both use the Cartesian basis 
or both use an angular-momentum basis for their momentum-space vectors, as is actually 
implicit in relation (§) between the Foldy and Lomont-Moses representations. We can then 
omit the factor T and replace the transformation flTjj) by 



Q F = D (0, 9, 0) O^D- 1 (0, 9, 0) . (19) 

The helicity representations of Lomont and Moses and of Shirokov are both seen to be 
unitarily equivalent to the standard Foldy representation, and the unitary transformations 
for both are rotations that rotate S3 into S-p. It seems surprising that the two helicity 
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representations appear so different. As Chakrabarti p4J points out, the Shirokov form of 
the operators J and K are of a singular nature that makes them appear unsatisfactory for 
many uses. However, the Lomont and Moses form (|T0| , |TT| ) also has first-order singularities 



at 9 = 7T, they are just better disguised. The differences in the unitary transformations ([5p 
and flT9| ) is clearer if U_ is expressed in terms of Euler angles. Since 

exp (— iS 3 (f>) S2 exp (iS^cp) = S2 cos — Si sin cf) = S • 0, (20) 

then from the Euler-angle form (|17]) we can equate 

U = D(M-<f>) , (21) 

that is, the unitary transformation U_ that rotates S3 directly to Sp is just the Euler-angle 
rotation D_{4>, 0, 0) required for the Shirokov form preceded by an axial rotation through the 
angle —(p. It is important to note that the transformations involved are spin rotations and 
therefore transform only spin components and the vector components of the photon states. 
We explore a generalization of such rotations further in the next section. 

While the position operators are not explicitly required to complete the Poincare algebra, 
they are implicit in the angular momentum and boost operators since, apart from internal 
degrees of freedom, it is the components of r x p and \{pv + rp) that generate rotations and 
boosts, respectively. One set of commonly accepted photon position operators consists of 
components of the Pryce operator p2]J18|, 



r p = ip a Vp~ a + (22) 

The parameter a in the expression 

p a Vp~ a = V + - [p, V] = V - ^ (23) 
p p 2 

depends on the integration weight used in the definition of the scalar product. Pryce used 
the parameter value a = 1/2, which is appropriate when the Hilbert-space states are electro- 
magnetic fields, and this value was implicitly assumed in Eqs.(0), (|TT1) , and (|l"3"|). However, 
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one needs a = — | when the vector potential |30| is used for the photon states. Alternatively, 
a definition of scalar product requiring a = can be used |T^J^1 . In the following, we leave 
a unspecified. 



The Pryce position operator rp is based on expressions in terms of a classical energy- 
momentum tensor of a noncovariant definition of the center of mass given by Fokker [|3T] . 
The components of rp are 3x3 matrices that do not commute with Sj. When expressed in 
terms of rp, the rotation and boost generators of Foldy [Eqs. (0) and ([?])] are partitioned 
differently into orbital and spin parts 

J = -pxr P +ppS (24) 
K = ^(r P p + pr P ). (25) 

While components of rp do satisfy the canonical commutation relations with p, they do 
not commute with each other Instead, they have a nonvanishing commutator that 

is analogous to the field of a Dirac monopole (but in momentum space) whose "charge" is 
given by the helicity operator pS fl"8|j . 

[r Pj ,r Pk ] = -ie jkl ^p-S. (26) 

As discussed in the next section, such monopole terms occur frequently in commutators with 
gauge potentials. 



III. PHOTON POSITION OPERATORS WITH COMMUTING COMPONENTS 

In I, a modification of the Pryce position operator was constructed to have commuting 
components . In this section we summarize the results of I, emphasizing how the new position 
operator can be expressed as a spin rotation of the weight-modified position operator for 
massive particles. We then extend the construction to an arbitrary functional dependence 
X P (0, (f>) among the three Euler angles 0, 9, x- 

The position operator 

r (0) = ip a V -a + _ a (0) p .g (2?) 
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with a(°) = 4>p^ cot 6 , was constructed in Ref. [|l9j by requiring that the components of 
have transverse eigenvectors in the directions 6 and <fi. The operator r(°) can be expressed 
as a spin rotation of ip a 'Vp~ ' : 

r (0) = iD (0, 9, 0) fVp^IT 1 (0, 9, 0) . (28) 

Since this has the same form as (|19D , it is equivalent to using ip a Vp~ a as the corresponding 
position operator in the helicity representation. The construction ( pq) may be understood 
by analyzing the problems of the Pryce position operator (22 ) . The reason components of rp 



do not commute is that the photon spin and its momentum are inextricably coupled through 
the restriction on allowed spin states, and in particular, the components of the momentum- 
space wave function depend on momentum direction. Consequently, the action of the usual 
momentum-space position operator for massive particles, zV, is no longer restricted to the 
argument of the wave function but also acts on its components. In expression (f2~8"|), on the 
other hand, the spin rotation D~ l rotates vector components of the state function from the 
lab to the reference frame, where p = e 3 , so that the operator iV can induce translations 
in its arguments without mixing its components. Finally, D rotates the translated vector 
back to the lab frame. This is a straightforward way to restore the usual translational role 
of z'V, and it is required to ensure that helicity eigenstates are not mixed by infinitesimal 
momentum-space translations generated by r^ -* away from the origin p = 0. This is an 
important property that follows from the invariance of the helicity subspaces under the 
Poincare group and has been confirmed, for example, in experiments demonstrating Berry's 



phase j34| for photons in curved optical fibers ||35|| . 

The term a^ ) = cj)p~ l cot 6 multiplying the helicity operator pS in the new position 
operator ( p7| ) is singular in the limits 9 — > 0, tt, that is, as p approaches the ±e 3 axes, and 
has the momentum-space form of the electromagnetic vector potential of a pair of Dirac 
monopole strings in coordinate space. The relationship is discussed in more detail at the 
end of Section V. It seems curious for the position operator to depend on the choice of 
coordinates in the laboratory. The singularity together with the dependence on laboratory 
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coordinates may have deterred other authors from including such a term. Nevertheless, we 
show below that some form of this term is required to give the correct phase changes of 
the photon state under rotations. Laboratory coordinates enter the formulation of r} ^ 1 in 
the spin rotation D_ from e3 to p and in the implicit dependence of p on its polar angles 
6, cf). The singular term arises from the differentiation of this rotation and from the path- 
dependent value of 0. It is consistent with the relatively large changes in that can result 
from infinitesimal changes in p at the string. In this section we explore a generalization of 
the rotation used in and show that the breakdown of the expected symmetry is related 
to the selection of a specific "gauge" for r . 

The spin rotation that plays a central role in this paper can be written in Cartesian form 
as the real unitary matrix of elements 

D jk = ej ■ e pk (29) 

that relate vector components in the orthonormal lab basis {e 1; e 2 , e 3 } to those in the photon- 
momentum basis {e p i, e P 2, e P 3} , where e P 3 = p. Thus, for any vector V, Vj = V • e 3 - = 
DjkV ■ e p k = DjkVpk ■ This gives the action of the matrix D as a passive transformation (V 
is fixed, the "observer" basis changes). However, only the relative orientation of the vector 
to the observer basis enters the formalism, and in particular, e^-e^ = Sjk = e P j-e p k. As a 
result, an active interpretation of the action of D_ (observer fixed, V changes) is also possible 
and often more natural. As an active transformation, D rotates each of the vectors ei, e 2 ,e 3 , 
of the lab frame into the corresponding vector e p i, e P 2,e P 3, of the photon p-frame, with both 
vectors expressed in the lab basis: e PJ = D_ &j . The mathematical consistency of the different 
interpretations is readily confirmed by writing out the column matrix representations of the 
vectors. Further details are given in Appendix B. 

The body rotation D_ generalizes the polar rotation used in I. One can express D_ in terms 
of spin matrices and the Euler angles 0, 9 and x as i n Eq. (|TT|), and it is given explicitly in 
Eq. (|B3|). The treatment in I corresponds to the choice x = 0, in which e p i = 6, e P 2 = tf>, 
and e p3 = p . For more general Xi 
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e p i = 0cosx + 4>sinx, (30) 
e p2 = 4* cos \~ Q sin X- (31) 

The angle % represents an axial rotation about p and is associated below with the phase 



of a photon with definite helicity. As seen by inspection of Eqs. (|30|J31|) , the transverse unit 
vectors e p i, e p2 , are just the polar unit vectors 0, <f>, rotated about p through the angle Xi 
so that use of the Euler angle x allows an arbitrary orientation of the two directions e p i, e p2 , 
in the transverse plane. The unit vectors describing states of definite helicity k, 

1 

e pK = -j= (e p i + ifte p2 ) (32) 
= e-^eg, (33) 

are just phase shifted relative to the helicity states e p °2 = (o + iK<pj / a/2 with x — 0- Under 
the general rotation (|TTD, the spin matrices themselves are transformed as 

Spi = DS 1 Dr l = SkD kj , (34) 

where in particular S p3 = S-p. Consequently, we can also express D in terms of rotations 
about the p-frame axes: 

D = DDD_- 1 (35) 

= e ~iD£^D-^ e -W£±D- 1 6 e ~iDhD- 1 x 

= exp (-iS p3 (p) exp (-iS p2 o) exp (-iS p3 x) ■ (36) 

Alternatively, we can make the axial rotation the last of the Euler-angle rotations by putting 
D = exp (— iSs(f)) exp (—iSzO) and noting 

D = D exp(-iS lX )D 1 D (37) 
= exp {-iDSzD^x) (38) 
= exp (— iS-px) exp (—iS&p) exp (— iS^O) . (39) 



The position operator that generalizes Eq. fl27|) is the 3x3 matrix 
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r = D (ip a Vp- a ) IT 1 (40) 
= ip a Vp- a -A, (41) 



where 



A = -%D (VZT 1 ) = i (VD) IT 1 (42) 

= ^P+ap.S, (43) 
p 

and a is a vector with dimensions of length. If \ is independent of 6, </>, then the axial rotation 
elements commute with V and r is given by as found in I. However, the momentum 
p is invariant under the axial rotation by x-, an d we can choose a functional dependence 
X = Xp (0> 0) • Then the matrix D_ is a function of the two parameters 9, <p, and direct 
calculation of the gradient of the rotation matrix gives 

a = a (0) + Vx P (44) 

with 

(0) = ^ = e 3 xpp-e 3 
P (e 3 xp) 2 

When operating on a subspace of definite helicity k = ±1, r is related to the Pryce operator 



r P of form (22) by a transverse displacement whose sign depends on k : 

r = r P — Ka . (46) 

This clearly manifests the basic coupling between the helicity (and therefore the spin) and 
the position. However, as we see in the following section, the eigenvalues of r do not depend 
explicitly on k . 

The matrix A, with transverse vectors as its matrix elements, has a simple interpretation. 
Writing 

VD = -i AD (47) 

and making the dependence of the rotation matrix D_ on p explicit, we find 
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D(p + dp) =H (p) + dp ■ VD (p) (48) 
= (l-idp-A)J2(p) (49) 

Thus, A dp generates the infinitesimal rotation induced when p is incremented by dp. 

The possibility of choosing different functions x p (6, 4>) dependent on the position 9, cj) in 
parameter space expresses the invariance of p under axial rotations as a local symmetry and 
is analogous to gauging a and hence r. The gauge group of axial rotations is U(l), and a is 
thus analogous to the electromagnetic vector potential. The scalar gauging function Xp is a 



concrete example of Berry's gauge [34| that sets the local phase of an eigenstate in parameter 



space. It also appears closely related to the function Go introduced by Staruszkiewicz 



in his earlier study [36] of parallel transport of photons. For all choices of single- valued 



differentiable functions Xp, the gauge field 

Vxa=-4, sinfl^O, (50) 

has the same monopole form. It follows that a itself cannot be replaced by a gradient term. 
As we discuss below (see especially Section V) however, x p is generally not single valued 
and can affect stringlike singularities in the field (|5*UD. The "abelian vector potential" a is 
just part of a nonabelian operator A = S x p/j) + ap-S . 

Berry |34| has found analogous monopole fields (|50|) both for the case of a pair of adiabatic 



states that are degenerate at an isolated point in parameter space and for spin-S* charges 
in a magnetic field that are degenerate where the field vanishes. Our derivation shows that 
such fields arise more generally in calculations of topological phase, even when there is no 
isolated point of degeneracy. Indeed, monopole terms of the form p/p 2 frequently arise in 
the commutators involving gauge potentials associated with massless representations of the 
Poincare group |I^j2^j32]j33|j37| -440| . In our case, direct calculation gives 



[ApAk] = ie m (V p x A) • e, (51) 
= ie jk i p-S 1 ■ e ; , (52) 

V p 2 v 



16 



which includes the monopole term p/p 2 multiplied by the helicity operator p • S. Similar 
terms also arise from the distinct "covariant derivatives" D' = V — ip x S/p 2 and D" = V — 
ka introduced for photon operators by Pati jllj and the Bialynicki-Birulas |p8|| , respectively. 
The common appearance of monopole terms and their associated string singularity (see 
Section V) is a direct consequence of the topology of momentum space [§§,0] (see also 
Appendix A). Here, the monopole "field" ( |50l) turns out to be just what is required to cancel 
the singularity in the commutator ( P6f) of the Pryce position operators. Furthermore, the 
relation (V + zA) D_ = implies that D = V + iA is the appropriate "covariant derivative" 
in our formulation. 

There is a wide range of possible gauge choices since Xp 0) can i n principle be any 
differentiable function. However, three choices of Xp are of particular interest: If Xp is the 
zero function, Xp = 0> then the transverse directions are the usual spherical polar unit 
vectors, a is just a^ ^ (|^), and r is given by Eq. (p7D . If Xp is Xp = ~ ( t ) -> then as noted in 
Section II, D_ is the Chakrabarti [51] transformation U_ and a takes the form 

a (i) = = pxe 3 

p sin p + p ■ e 3 

in Eq.(^TJ). We note that this choice eliminates the singularity in a at = 0, but doubles 

(2) 

the strength of the one at = n (see also Section V). If instead Xp is Xp = —<ficos0, then 
a is given by 

a ( 2 ) = -0sin#. (54) 
p 

In this case, both string singularities are removed, but a is now path dependent, since the 
value of depends on how many times the path wraps around the e3 axis. Consequently, 
a ( 2 ) j s "nonintegrable" ||38|| . By subtracting another path-dependent term, it is also possible 
to reintroduce a string singularity about a different axis. 
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IV. POINCARE GENERATORS AND COMMUTATION RELATIONS 



Here we examine the position operator r in the context of the Poincare algebra and 
write explicit momentum-space expressions for the localized bases. To see the relationship 
of our position operator r and the Poincare operators in the existing literature, note that 
translation generators p and H = pc are unchanged, and the Foldy J and K operators can 
be written in terms of an orbital angular momentum operator 

L (r) Erxp (55) 

corresponding to our position operator (^TJ), and an effective spin operator, which has a 
transverse contribution in addition to the helicity term that is frequently assumed to be the 
total spin contribution 



S (r) E(axp + p)p-S. (56) 



We obtain 



J = L (r) + S (r) (57) 



and 



K = ^ (pv+rp) + p x S (r) (58) 
= \ (pi+IP) + ap • S (59) 

as the Foldy operators (generalized to arbitrary a ) in terms of our position operator. With 
these expressions, transformation to one of the helicity forms is straightforward. By inverting 
Eqs. (|40|) and fl34|) the position operator becomes i2 _1 ri2 — ip a ^7p~ a , while i2 _1 P ■ SID =S^. 
The general helicity representation is given by O h = D^OD, which gives 

J h = (tp a Vp- a ) x p+ (a x p + p) 5a (60) 
K h = % -p a (pV + Vp) p~ a + SipSs (61) 
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If Xp = 0) then a = (f>cot9/p and Jh and Kh reduce to the Shirokov operators given 
by Eqs. (|12"D and C [3~3|) generalized to arbitrary a, while if Xp = ~ ( t ) i a is given by 
cf)(cos6 — 1) / (psin#) = — </>sin#/ (1 + cos#) and Jh and Kh reduce to the Lomont and 
Moses operators, Eqs.([T0|) and ( pTTj) with arbitrary a. Thus both the Lomont and Moses and 
the Shirokov transformations are special cases of Euler angle rotations relative to the Foldy 
form. In any irreducible helicity representation the position operator is just ip Q Vp _a as it 
is for massive particles in the Foldy form, since the transverse directions are fixed in the 
reference frame where the photon momentum is parallel to e 3 . The helicity reference frame 
is the photon frame, while the Foldy frame is the lab frame. 

The commutation relations satisfied by the generators of time and space translations, 
rotations and boosts, H = pc, pt, Jj and Kj, are the standard ones required by the Poincare 



algebra: 

Kj,Pk 



Jj i J k 

■ i5 jk H/c, 



K 3 ,H 



Kj,K k 



Jj,Pk 



ICjklPl, 



icpj, 



J j: H 



0, [pj,H] = 0. Analogous commutation 



relations are valid if is substituted for Jj and (prj + rjp^j /2 for Kj, since they are 
unitarily equivalent to the zero-spin case of the above. Commutation relations involving 
Tj in the lab frame can be most simply derived by performing a matrix rotation to the 
photon reference frame, in which r is replaced by ip a Vp~ a and noting that components pj 



in the momentum representation commute with the rotation matrix D_. Thus 



rj,n 



rjiPk 



i5 



Tj,H 



iHpj/p 2 



3 1 1k 



Tj , 'S'pfc 



and 



= 0, 
0. It 



follows from the last two relations that a photon can simultaneously have definite helicity and 
either a definite position or a well-defined spatial-component of orbital angular momentum. 
In other words, measurement of either the position or the orbital angular momentum of a 
photon state does not change its helicity. Furthermore, because 



fep-S] =0, [K fc , P -S] = 0, [ Pfc ,p-S]=0, 



(62) 



the helicity is an invariant of the Poincare group; it is invariant under all rotations, boosts, 
and translations (as long as these avoid the origin p = 0). Consequently, every representation 
of the Poincare group for photons can be reduced to the direct sum of representations for the 
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two helicities, and every irreducible representation will be carried by states of a single helicity. 
It may therefore be convenient to specify the operators for the separate invariant helicity 
subspaces. The momentum and energy operators are unchanged, the position operator is as 
given above in Eq. (^6|), and J and K take exactly their Pryce forms 



J = L (r) + S (r) = r P x p + f>K 



(63) 



and 



K = ~ (pr +rp) + a k = ^ (pr P + r P p) 



(64) 



with rp = r + (ta . 

It is important that our position operators obey the correct dynamical equations. In the 
Heisenberg picture, the dynamics are determined by the equation of motion, which from 
Eq. (gDp is 

dr dr 



dt dt 



+ i [H, r] = i[cp, r] = cp . 



(65) 



Thus our theory predicts that the photon has a velocity v = cp, as required. We also note 
that r (f41~D is Hermitian and symmetric under time reversal. One may in addition expect r 
to change sign under parity inversion, but this depends on the gauge potential. It is valid 
for r^ ), but for other gauges we must generally replace the inversion in x p by the gauge 
change Xp ~~ Xp i n order to ensure the invariance of Vx P • 



In general, the position operator rotates as a simple vector under L/ r ) since 
itjkili , but not under J : 



3 i 



Sj , T_k_ 



(66) 



From 



P-S, r 3 



and Eq . (|56|) , we thus find 



itjkiri + [(a x p + p) • ej,rk\ k 

d , 



(67) 
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The extra term represents a deviation from the usual commutator (0), and is due to the 
coupling of the momentum and spin of a photon. Note that the deviation vanishes for 
rotations about p. It also vanishes in the special case of a rotation about e3 when Xp is 
independent of (f>. 

For a photon with helicity k the momentum-space states localized at r are 

%, iK (p) = Np a e- ir ' p e pK , (68) 



as may be verified by direct application of the operator r fl4Up. In particular, the localized 
state at the origin, 

*o, K (p) = Np a e pK (69) 



gives 

L (r) ep« = -p x re^ = (70) 

since ^o,k is an eigenvector of r with eigenvalue 0. The relation r e pK = also follows directly 
from the transformations (f40"|) and ( |B1| ). Using J= + S ( - r ' ) where is given by Eq.([5^), 
we find the rotation by an infinitesimal angle <i£ about the axis £ to be 

exp (-zJ • d£) ep* = exp (-iS (r) • d^j (71) 
= exp (— m(a x p + p) • d£) e pK 

with d£ = £<i£ . Thus the rotated transverse basis vectors of helicity n satisfy 

exp (-iJ • dg) e pK = exp [-m (d\ - dx P )} e^, (72) 

and are thus changed only by an infinitesimal phase shift 

— k (dx — dxp) = — ft(a x p+p) • d£. (73) 

Here, dx P = dp ■ V% p is the change in the function % p (9, (j>) that results from the change 
dp in the photon momentum. (Note that a and depend on the orientation of p, and 
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therefore ordered integral expressions of the Dyson type are required for finite rotations. See 
also the Appendix.) 

Since e p0 = p is a longitudinal basis vector corresponding to helicity k — 0, it follows 
from p- S e p0 = that Eq. ([7^) is also true of the longitudinal vector. Indeed, Eq. ( |72| ) may 
be considered an extension of the expected invariance of the radial vector field p under 
rotations. Also, if £ = e 3 , and Xp is constant, then it is readily shown that d\ — d\v = so 
that this cylindrically symmetric transverse vector field is invariant under rotations about 
e3. Note further from Eq. (|49|) that with dp = d$, x p, the infinitesimal rotation factor 
becomes 

H (P + dp) IT 1 (p) = 1 - iA ■ dp = 1 - i (S - S (r) ) • d£ . (74) 

The part — iS ■ d£ gives the additional rotation whereas iS^ ■ d£ corrects the axial rotation 
implied by the functional dependence Xp ($> 0) • 

One can understand the extra term in the commutator of the position with the angular 
momentum by comparing Eqs.(|67D and (0). Rotations generally change the "gauge po- 
tential" a in a way analogous to transformations of the electromagnetic potential A under 
Lorentz transformations. The <9(a x p + p) • ej/dpi term in Eq.(p7|) is required in order to 
give the axial spin rotation needed by the matrices associated with each component of r. 
Such axial rotations result in the phase change of eigenvectors during a rotation that main- 
tains the functional dependence of Xp an d is required in order to give the correct Berry's 
phase for photons (see the next section) . The possible transformations are discussed in more 
detail in the Appendices, and Eqs.([73D and (|67D are independently derived in Appendices B 
and C respectively. 

V. RELATION TO BERRY'S PHASE 

The phase angle — k (dx — dx P ) is important and has measurable consequences. In par- 
ticular, it can be integrated to give the total phase change of a state vector of helicity k 
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when transported along a closed loop in parameter space. The result depends on the type 
of transport, but it will not depend on any single- valued functional choice of Xp- However, 
it is natural to make Xp multi- valued and path-dependent (nonintegrable) , and in that case 
the phase change will depend linearly on the change A% p around the closed loop, as we 
show explicitly below. If we use parallel transport, the angle — k (dx — dx P ) is Berry's phase 
3~4] , |43H j K , a topological phase accumulated by the photon. This is consistent with a more 



general association of Berry's phase in quantum systems with geometrical angles. p4| 

While Berry's derivation assumed adiabatic transport of energetically discrete states, his 
results can also be applied to degenerate helicity states in parallel transport |4^,|35| . Parallel 
transport on the spherical surface of constant radius p is most easily realized by piecing 
the path together from many small segments of great circles and employing nonrotating 
transport along each segment. Each great-circle segment requires an axis of rotation that 
is perpendicular to p but generally changes as p moves. From Eqs. (|73|) and ([72]), the 
accumulated phase for great-circle segment (p ■ d£ = 0) is 



- k (Ax - Ax P ) = -k J a x p ■ d£ 
= —k <j) a • p x d£ 

= k (j> a-dp , (75) 
where we noted dp = d$, X p . The same result (|75|) can be obtained directly from Berry's 



derivation [[34| of the geometric phase j K . We can write his starting point in differential form 

as 

dl K = i(u K (p) |Vu«(p)) dp (76) 

since the relevant parameter for the photon as it is guided in an optical fiber is its momentum 
p . In our case, the eigenstate u K (p) is proportional to the column vector 

ep« = D (77) 

so that Berry's relation ( j76|) reduces to 
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d lK = i (ejD} (V-D) e«) • dp = (^^A^) • dp (78) 



pK f A e pK ) ■ dp = /ta • dp (79) 



where the matrix sandwich indicated by angular brackets (• • • ) is a vector- valued lxl matrix, 
and we have noted that Sxp has a vanishing diagonal in the e pK basis. In terms of the polar 
angles 9, <p, displacement on the spherical surface is given by dp = </> p sin 9 dip + 0pd9 , and 
therefore with Eqs. ( f44| , f45|) for a , the line integral (|75| ) gives 

- K( A x -A Xp )= K (/cos0# + /v Xp . (i p) 

= k (<j) cos9d<p + Axp^j . (80) 

The result depends on the function Xp • We must choose a function that eliminates 
superfluous axial rotations. Recall that the basis vectors e pK are obtained from the reference 
frame by the rotation D_{4>, 9, x P )e K . To avoid unwanted axial rotations, we ensure that p is 
also parallel transported by this rotation. This requires the function choice 

x P (M) = x P 1} (M) = -</> (81) 

to make D_{4>, 9, Xp) equivalent to the direct great-circle rotation from e3 to p. The phase 
shift is then 



7 K = — k (A X - Ax p 1} ) — cos 9 d(p — lix^j = -kM, (82) 



where Q is the solid angle enclosed by the loop. This is Berry's geometrical phase |34| which 
has been confirmed in experiments on light in helically wound optical fibers [p5f| . The sign 
convention |34| is taken such that the dynamic phase of a stationary state (due to a factor 
exp (— iuit)) decreases in time (see also Appendix D). 

One often seeks to generalize the above result by an application of Stokes theorem. In 



our case 



^"a (1) -dp = J (V x a (1) ) ■ pp 2 dQ 



- / p ■ ^p 2 dil = -Q , (83) 
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where f2 is the solid angle of the integrated area. The result should be valid as long as the 
integrated area and its boundary avoid singularities and branch cuts. Since it is the curl 
of a that appears in the surface integration fl83|), gauge transformations a — > a + Vxp will 
not change the result if Xp is a single-valued function of p. This is consistent with the line 
integral fl8HD since for any single- valued Xpi the difference A% p vanishes over a closed loop. 
If, as in the optical- fiber experiments [ 35f] , p describes a circular path that makes a fixed 



angle 9 with any fixed direction, then Q = 2tt (1 — cos 8) so that the change in phase angle 



is 



1k = K j> dp ■ a (1) = 2ttk(cos# - 1) . (84) 

However, the derivation also apparently works if is replaced by a^, even though these 
results should differ by kAxp = ~~ ^ 7CK f° r a closed loop around e 3 . While such a difference is 
not observable, we can trace its origin to the nonintegrable nature of xp an d the associated 
string singularity through the integrated area of the surface integral (^3|). The argument 
can in fact be turned around to imply that because the difference cannot be observable, the 
helicity must be quantized. 

Since the string contribution is frequently ignored, it may be useful to elaborate its role. 
Recall that Dirac's magnetic monopole contains a string to bring the magnetic flux to the 
monopole inside a solenoid with a vanishing diameter. Although strings are not evident in 
our expression fl5"0|) of V x a, their presence is implied by the expression (|44|) for aS ^ . The 
physics is clearer if we write a^ ) as the limit of a nonsingular function: 

a (0) = 4> lim ^ ? , 85 

where pj_ = psin^ is the momentum-space distance from the e 3 axis and po represents the 
approximate diameter of the solenoid. The "magnetic field" corresponding to the "vector 
potential" is 



V x a (0) = lim 



sin^ 9 2p 2 Q cos 9 



-P _2 , _2 + e 3 



(86) 
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The first term becomes the monopole field in the limit, and the second term represents 
the two axial strings. It vanishes everywhere except on the axis and can be expressed as a 
two-dimensional Dirac delta function, giving 

V x a<°> = + 27re 3 cos 95 {2) (pj , (87) 
pi 

where we can also write 5^ (p±) = S (pi) S Q02) , and the coefficient of the second term on the 
rhs has been chosen to give the correct surface integral over p±dp±d<p at fixed (pcos^l ^> po- 
On the string, cos 9 = ±1, so that we have two half strings along the two halves of the 
e3 axis, both taking "magnetic flux" away from the origin. The string term adds exactly 
27r to the surface integral (0), thereby restoring Stokes theorem and bringing the surface 
and line integrals into agreement. As mentioned above, gauge transformations with various 
functional forms Xp 0) can shift the strings and replace them by explicitly nonintegrable 
functions. The choice of Xp makes a nonsingular over the integrated surface in Eq. fl33|) 
between the loop and the upper pole 9 = 0, but a^ ) has a singularity penetrating the same 
surface, giving both line and surface integrals for the phase difference — k (Ax — Ax p ) that 
differ by kAx^ ■ The agreement also extends to use of the surface that includes the pole at 

(2) 

9 = 7r for both cases and 1. For the choice xp the line and surface integrals both vanish if 
the branch cut is inserted explicitly in the surface integral, so that the nonintegrable function 
in Eq. (|54] ) is replaced by 

a< 2 ) = -[</>- 2nh (0 - O )1 sin 9 , (88) 
p 

for example, with < 0o < 27r where h((j) — 0o) is the Heaviside step function. 

VI. DISCUSSION 

In this Section, the nonexistence proofs and recent papers concerning localized states 
are briefly reviewed, the new photon position operators are discussed in the context of this 
literature, and our results are summarized. 
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The most quoted paper is that of Newton and Wigner [fL3|| . These authors assumed 
a rotationally invariant set of localized states and arrived at the position operator of the 
form ip a Vp _a for spinless particles with or without mass. They also obtained an expression 
of the position operator for massless particles of spin |. Regarding photons, they stated 
that for S = 1 and higher "we found that no localized states in the above sense exist. 
This is an unsatisfactory, if not unexpected, feature of our work." As a result of their 
conclusion, it is frequently stated that a (spatial) photon wave function does not exist [||. 
The localization postulate adopted by Newton and Wigner is strong (any displacement of a 
localized state is assumed to make it orthogonal to states of the undisplaced set) and has 
been the focus of a number of more recent studies. In particular, Wightman and others 



14| have used generalized imprimitivities ||5| to reformulate localization more rigorously in 



terms of localizability in a region. However, their work did not alter the conclusion that 
a single photon is not localizable. Other authors have sought effective wave functions that 
satisfy a somewhat relaxed localization condition P JT^j27|j4^j47j1 . 



We suggest a different potential problem with the conclusion of Newton and Wigner 
for massless particles of spin S > 1. To ensure a rotationally invariant linear manifold 
of localized states for a system with total angular momentum quantum number j , they 
assumed a complete set of 2j + 1 wave functions ipjm, — 3 < Tn < j '• , where m is a component 
referenced to an external direction. While the existence of a complete set is sufficient to give 
a rotationally invariant manifold, it is not necessary for massless particles of spin S > |. 
Massless particles with spin have only two spin states, namely those corresponding to the 
helicities ±S. For a system of states at the coordinate origin, the orbital angular momentum 
vanishes and j = S. The states in the linear manifold are characterized by components of 
j not along a space- fixed direction but along the momentum direction p. For S > |, the 
manifold is not complete and consequently it cannot describe a state with spin quantized 
along an arbitrary direction. However, it can describe the allowed states with either helicity. 
Furthermore, since the helicity operator commutes with the generator J of rotations, the 
two helicity subspaces are separately rotationally invariant. Because the helicity eigenstates 
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form a complete rotational set only for S < ~, it is clear why the Newton's and Wigner's 
insistence on a complete rotational manifold is stronger than necessary for massless particles 
with 5* > 1 . 



In the papers based on the method of generalized imprimitivities it appears that 
the system of commuting imprimitivities, through which the position operator is defined, 



is assumed to be independent of spin. (This arises from Wightman's axiom V [|Tj]). As 
seen above, for example in Eq. (|46|) , the spin is inextricably coupled to the momentum and 
thereby to the position operator. Position operators are known to exist for massless particles 
of spin 0, and it may be possible to find a system of imprimitivities like that for the scalar 
case for each value of the invariant helicity. In any case, the imprimitivities must be distinct 
for different helicities because the position operators are. 

In several more recent proofs of the nonexistence of a photon position operator with 
commuting components, operator algebra was used with the assumption that the position 
operators satisfy Eq.(|J) |I5| - |I8fl . However, we have shown above and in Appendix C that 
the matrix r includes a "gauge potential" a that is transformed by rotations. As a result, 
the commutator [Jj,rk] contains an additional term involving the spin, and indeed that 
this extra term is required in order to give the correct phase of the rotated photon state. 
The components of the position operators ( f40|) thus satisfy Eq.(^), which unlike the more 
familiar relation (|3]), is compatible with commuting components. Thus, nonexistence proofs 
that assume Eq. (§) do not apply. 

Our construction of the photon position operator uses a spin rotation to decouple the 
spin from the momentum while the gradient operator acts to generate momentum-space 
translations. It seems to be the most natural way to retain the usual role of the position 
operator in momentum space. While we have not touched on important questions about the 
limits to which the photon position is observable in a dynamic measurement, our demonstra- 
tion that a photon position operator does in fact exist means that there is no nonvanishing 
commutator [r^rjj to limit the calculation of photon probability amplitudes. Simultane- 
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ous eigenvectors of the operators are available for calculation of the probability that the 
corresponding eigenvalue is observed. Of course limitations arising from Fourier analysis, 
similar to those applicable to massive particles, still apply. Our work therefore supports 
the view that photons have wave functions that are not qualitatively different from those 
of massive particles, as concluded by Bialynicki-Birula jRJ and Sipe ||. The rules of quan- 
tum mechanics require that each observable be represented by a Hermitian operator. The 
localized basis sets found here makes it possible to treat photons like massive particles in 
quantum calculations of interference experiments and other situations where particle am- 
plitudes in coordinate space are useful; in both cases the usual rules of quantum mechanics 
can be applied. 

There is no unique representation of the photon position operator. Just as there are many 
spin bases that can be used to describe the internal state of a massive particle, there are 
many bases that can be used to describe the combined translational and internal motion of a 
photon. The peculiar aspect of the photon is that, because of the coupling of the spin to the 
momentum, the position operator is a matrix that does not commute with the spin. Different 
position operators with commuting components can be defined by specifying the dependence 
X P {0, 4>) of the axial rotation angle on the polar angles of p. Different choices of Xp l ea d 
to different "gauge potentials" a for the phase of the photon state at different positions in 
momentum space. In general, phase changes depend on the path and are thus described 
by a nonintegrable function. While a specific choice of gauge usually results in operator 
expressions with less than the expected symmetry of the Hamiltonian, the symmetry is 
obtained within the full group of possible gauges. The forms that result from a selection 
of Xp are ah unitarily equivalent to each other and to the position operator for a massive 
particle, zp a Vp~°, where the unitary operator is the rotation through Euler angles. The 
operator ip a Vp~ a has eigenvectors with a fixed direction in space and is independent of 
the spin. Consequently, for massive particles the description of spin is a separable problem. 
However, the spin and momentum of a photon are inexorably coupled, since the direction of p 
determines the direction of the observable component pS of the internal angular momentum. 
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Position operators with transverse and longitudinal eigenvectors could in principle be used 
to describe a massive particle, but this is probably not useful since ip a Vp~ a is a simpler 
alternative. For a massless particle this choice does not exist, reflecting the fact that for a 
photon, the orbital and spin angular momenta are not separable. 

In summary, the arguments presented here show that a photon is much like any other 
particle in that its position is an observable described by a set of three commuting Hermitian 
operators. However, the photon (as well as other massless particles of spin S > |) has 
only two linearly independent spin states, and in these states the spin is coupled to the 
momentum. As a result, its position operator is a matrix that does not commute with the 
spin. Different selections of the function x P (9, 0) generally give different position operators, 
so that the position operator is not unique and does not transform under J as a simple 
vector. However, the eigenvectors of any one of these unitarily equivalent position operators 
gives a basis of localized states with unique eigenvalues that are independent of helicity, and 
there is consequently no disagreement as to the actual position of the photon. Contrary to 
the traditional view, localized basis sets do exist, and it appears that photon wave functions 
can be defined according the usual rules of quantum mechanics. 
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APPENDIX A: DECOMPOSITION OF ROTATIONS 



The purpose of this appendix is to derive relations between an infinitesimal rotation about 
an arbitrary axis in three-dimensional space and angular parameters in a product of rotations 
around specified axes. These relations are then used to predict rotational properties of the 
photon position operator and its eigenstates. The angular parameters are the Euler angles 
that specify the polar and azimuthal coordinates 9 and 0, respectively, of a given direction p 
together with an axial angle \ about p. One can express the Euler-angle parametrization in 
terms of active rotations about spaced-fixed axes: with the p-frame axes initially coincident 
with the space-fixed axes e 1; e 2 , e 3 , of the lab frame, a first rotation by x about e 3 is followed 
by a rotation by 9 about e 2 and finally a rotation by about e 3 . The axial rotation angle 
X does not affect the direction p and may be chosen to be a function, say Xp ($> 0) > °f the 
"local coordinates" 9, that specify p. A rotation by the Euler angles (0, 9, Xp (#>0)) has 
only two degrees of freedom and is uniquely determined by the direction p. 

The derivation is easily constructed using tools of Clifford's geometric algebra C£ 3 of 3- 
dimensional space |B§-[K|. Those not familiar with the algebra can follow the derivation in 
terms of the 2x2 matrices that form the common matrix representation in which the Carte- 
sian unit vectors ei, e 2 , e 3 , are replaced by the corresponding Pauli spin matrices familiar to 
physicists. In the algebra, vectors are rotated by transformations of the form 

v -> RvR\ (Al) 

where R (£) = exp (— i£/2) G SU (2) is the element for a rotation by the angle £ = |£| about 
the axis £ . The rotation specified by the Euler angles (0, 9, x) is given by the rotation 
element 

11 (0, 9,x) = R (0e 3 ) R (0e 2 ) R ( X e 3 ) (A2) 
= exp (-z0e 3 /2) exp (-i9e 2 /2) exp (-zxe 3 /2) . 

An additional rotation will generally change all three angular parameters 0, 9, x ■ We want 
to determine the changes caused by an infinitesimal rotation R (d£) . 
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To clarify our objective, we first consider the simple case in which d£=d£ e 3 . The only 
effect of R (d£) is to increment the azimuthal angle '■ 

K (0, 9, X )^R (d£e 3 ) K (0, 9, x) = K (<P + 9, X ) . (A3) 

This result is easily written in terms of the rotations 71 (0, 6, Xp (9, 0)) that have only two 
degrees of freedom: 

R(d& 3 )K(4>,6,x P (e,4>)) (A4) 

= n (0 + d£, e, xp {o, + dO) R (-e 3 ^ P ) , 

where x P (0, + d£) = x P (9, 0) + dx P and dx P = d£ dx P / d<fi . 

Now we generalize this approach to an arbitrary infinitesimal rotation by d£ = d^ei + 
d&e2+d£ 3 e 3 .We initially consider x an independent parameter and solve R (d£) 1Z (0, 9, X ) = 
1Z (0 + d(j), 9 + dd, x + dx) for the infinitesimal changes d<p, d9, dx in the Euler angles. To 
first order in the changes, !Z((f> + d(p, 9 + d9,x + dx) 



1 - ~ (e 3 d(p + e- ie3,l, e 2 d9 + Re 3 R)d X ) 



n{<j>,9,x) 



1 

1 — — [ei (sin 9 cos <pdx — sin <pd9) + e 2 (cos <pd9 (A5) 

+ sin 9 sin 0c?x) + e 3 (d<f) + cos 9dx)]} ft (0, 9, x) , 

in which we have explicitly accounted for the lack of commutivity of rotations about different 
axes. Equation ( |A4j) therefore implies 

d£ = ei (sin 9 cos (f>dx — sin (pd9) (A6) 
+e 2 (cos 4>d9 + sin 9 sin 0g?x) + e 3 (dtp + cos 6*c/x) , 



which is readily solved to give 



d^ 3 — cot 9 cos + g?£ 2 sin 1 



d9 = d^2 cos — d£i sin (A7) 

, d£i cos + <i£ 2 sin0 

OX = — a • 

sin 
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Note the singular nature of the relations for d(j) and dx in the limit 9 — > . Singularities are 
common whenever general rotations are parametrized in terms of rotations about specified 
axes. They are related to the nonuniqueness of the parametrization for some rotations, for 
example all of the rotations D (<p,0,x — 4>) for fixed \ and arbitrary <fi are equal. We can 
also express the rotation (A5) in terms of rotations with two degrees of freedom 



R (d£) U (0, 9, Xp ) =n(<P + d0,9 + d6, Xp + d x ) 
= n(<i) + d<j),8 + d9, Xp + dxp) R ({dx - dx P ) e 3 ) 
= R ((dx - dxp) p) H (0 + d(j), 9 + d6 } Xp + dx P ) 



where x P = X P 



^) and 

dx - dxp = dx - 
cos « 



dx 



+ 



06 



+ 



sin 9 
sin <p 
sin^ 

dCs 



1 + cos V^f- I + sin <P-q^- 

dXp 



1 + cos 9 



dXp 



COS) 



09 



dCi 
d& 



dXp 



(axp + p)'^, 



(A8) 



(A9) 



(A10) 
(All) 



where in the last step we used the definition of a . While the derivation employs tools 

of geometric algebra, the result is generally valid for spatial rotations in three-dimensional 
space. 

A similar analysis gives the result of an infinitesimal rotations added to the rhs: 



K (</», 9, X p) R (dg) = TZ(<P + d<p,9 + dd, X p + d X ) 
= n(<p + d<f),9 + d9, Xp + dxp) R((dx- dx P ) e 3 ) 



(A12) 
(A13) 
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In this case, we obtain 



-d£[ cos x + d£ 2 sin x 



sin 9 

dO = d£[ sin x + dg 2 cos x 
dx = d£' 3 — dtp cos 9 
dx ~ dx P = d£[ 



(A14) 



cosx ( a . dx P 
cos 9 + — — 



sin^ 



sin 9 



cos 6> + 



smx 



cosx 



dx P 
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dx P 
89 



(A15) 



These results can also be obtained from the Hermitian conjugate of the relations (A5) and 
([A3]) by noting that (0, 9,x) = D (— x, —9, — 0)- Such added rotations from the rhs are 
rotations in the p-frame and can also be expressed by 



K (0, 9, xp) R {dtfe) = R (d^e pj ) U (0, 6, Xl 



(A16) 



Evidently the Cartesian components of d£' are just the p-frame components of d£, that is, 
d& = d£ P j . We make use of this result in Appendix D. 

APPENDIX B: APPLICATION TO ROTATIONS OF STATE VECTORS 



The helicity state vectors e pK are obtained by rotating the constant column vectors e„ 



by the rotation matrix 



e pK = J2((f),9,x)e H 



D = exp (— iS 3 (f) exp (—iS 2 9) exp (—iS 3 x) 



generated by the spin-one matrices Sj. Starting in the Cartesian basis with 



e 







, e ± i 



V2 



±i 

v° / 



(Bl) 



:b2) 



(B3) 
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the rotation matrix is explicitly 
/ 



D 



cos 9 cos cos x — sin sin x — (sin cos x + cos 9 cos sin x) sin 6* cos 
cos 6* sin cos x + cos sin x cos cos x — cos 6* sin sin x sin 6 1 sin 



\ 



V 

and we thus find 



sin 9 cos x 



sin sin x 



cos 9 



J 



(B4) 



sm9 cos0 



sin^ sin i 
cos 6* 
/ 



(B5) 



/ 



-IKX 



V2 



\ 



\ 



« = ±1. 



(B6) 



/ 



cos 6* cos (f) — ik sin i 
cos 6* sin + in cos 
— sin# 

Note that e p0 is a matrix representation of p, and e pK depends on x simply through the 
phase factor e~ tKX . We put x — Xp 4>) to obtain vectors e pK that depend only on the two 
angles 9 and 0, and hence only on the direction of p . The direction of p can be changed by 
a further rotation of e pK by d£ of the form exp (— id£-S) e pft 



D((f) + d<f),9 + d9, x P + dx)e. 

exp [-ik (dx - dx P )] D((f) + d(f),9 + d9, x P + dx P ) e« 
exp [-ik (dx - dx P )} e p+dpK . 



(B7) 



(B8) 



The rotation matrix exp (— id£-S_) shuffles the components of the vector e pK . A complete ro- 
tation of e pK also includes the operator exp (—id£ • L) which changes the angular arguments 
9, so as to rotate p to p — dp . We thus find exp (— id£-J_) e pK 



exp [-m (dx - dx P )] exp (-id£ • L) e p+dpK 
exp [-in (dx - dx P )] e pK 



(B9) 
(BIO) 



and consequently the only effect of a rotation of e pK by dC, is to change its phase by 
— k (dx — dXp) , where dx — dx P is given by Eq. (A10). This agrees with the partition 
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of J into 



since 



and 



J = L (r) + S (r) (Bll) 

L (r) ep, = D LD~ V = D Le« = (B12) 

S^=p-S(axp + p) (B13) 

p-S ep« = DS^D' 1 ^ = DSz = Ke^ (B14) 

(a x p + p) • d£ = dx — dx P • (B15) 

APPENDIX C: APPLICATION TO ROTATIONS OF THE POSITION 

OPERATOR 



The new position operator has components 

d 

n = ip a D—D_- 1 p- a . (CI) 
- dpi 

The components of vector operators are expected to rotate into one another as given by the 
commutation relation (^). However, an additional factor arises by the need to transform 
the matrix associated with each vector component by an axial rotation. The required spin 
rotation is given by that for the rotation matrix D_. From relation (A8) we see that the 
rotation transforms D (0, 9, % p ) to (1 — iS ■ d£) D ((f), 9, x P ) 

= D{<p + d<f),9 + d9, Xp + d Xp ) (1 - 1S3 {d X - d Xp )) (C2) 
= (1 - ipS (d x - dx P )) D(4> + d<f),9 + d9, Xp + d Xp ) , (C3) 

where from Eqs. (All) and (A10), 

P-S (d X - d Xp ) = P-S (p + a x p) = S (r) • d£ . (C4) 

The result is to add S/ r \rk to the rhs of Eq. (^. The added term is required for 
consistency with the phase changes (BIO) induced in the basis vectors by rotations. 
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APPENDIX D: ROTATIONS OF THE REFERENCE FRAME AND BERRY'S 

PHASE 



An extra rotation acting on the rhs of R (0, 9, x) is equivalent to the opposite rotation to 
the reference frame. As shown in relation ( |A16| ), it is also equivalent to a rotation from the 
left by the same angles but about p-frame axes instead of lab axes. This relation is useful 



in deriving an expression for Berry's phase [B3]. Berry's phase is a topological phase that 



arises from adiabatic transport of a discrete state around a closed loop, and it is usually 



derived as a purely quantum phenomenon . For the polarized light, however, it 



was shown by Haldane |51] to be a consequence of classical differential geometry. We derive 
it here by means of the classical rotation operators introduced above. 

Transport in the adiabatic limit corresponds to parallel transport in the given gauge 
|35| , |52[| , which is equivalent to nonrotating (Fermi- Walker) transport on geodesies. We 
achieve this in rotations of p by piecing together the closed loop from a sequence of great- 
circle rotations, in each of which the rotation axis is orthogonal to p and the axial rotation 
vanishes. Thus, in Eq. (A8) we vary £' such that d£' 3 = . Then from Eq. (A14), 

dx = — cos 9d(fi (Dl) 

We vary the parameters and £' 2 so as to complete a closed loop, increasing by 2n. The 
condition d£ 3 = ensures that the change in p is always along a great circle. There is a 
change in the axial angle — x by 

- Ax = I cos 9 d(p (D2) 



even after p has looped back to its original direction. If the initial orientation of the photon 
frame is given by R ((f), 9, x) with the momentum direction p, the orientation after has 
increased by 2n is R (0 + 2n, 9,x + A%) , giving a the net rotation of 

i?(0 + 27r,fl, X + A X ) iT 1 (0,0, X ) 
= exp [-i (x + Ax) p/2] exp (-i7re 3 ) exp (ixp/2) 
= expH(A X + 27r)p/2], (D3) 
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where we used the symmetry 

R (0, 9, x) = exp Hxp/2) R (0, 9, 0) , (D4) 

which allows us to cancel the 9 rotations, and we noted that exp (— iires) = — 1 = 
exp (-itp) . The total rotation under parallel transport about the loop is thus the rota- 
tion about p by the angle rj = Ax + 2n = § (1 — cos 9) dep. 

Note that the result does depend on the parallel transport of p. If we rotate p directly 
around the space-fixed e3 axis, then we simply change by 2n and the axial angle \ does 
not change. The mathematical confirmation of this result is found from Eq. (A7) with 
d£i = d^2 = 0. The rotation in this case is not on a great circle (unless cos 9 = 0) and p is 
therefore not parallel transported. 

The phase change for electromagnetic plane waves of given helicity is proportional to the 
rotation angle. To see the relation, consider electromagnetic plane waves of helicity k, for 



which the field F = E + icB is 49 



F = (1 + p) E exp [in {ut - p ■ x) p] . (D5) 

An additional axial rotation by rj transforms this to 

exp (— irjp/2) Fexp (ir/p/2) = Fexp (irjp) (D6) 
= (1 + p) E exp [in (out — p • x + Krf) p] (D7) 

and gives the phase change — kt] relative to —cut. This is exactly Berry's phase 



k (J) (cos 9 — 1) def) = — nVl , (D8) 

where fl is the solid angle subtended by the loop. It is just that calculated directly from 
the line integral in Section V when Xp — Xp = ~ ( P so that the rotation matrix D (0, 9, x P ) 
rotates S 3 directly into S • p, and it agrees with the relation between angle and geometric 
phase found by Berry [ |53|| . 
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